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Abstract 

Multiple orthogonal polynomials (MOP) are a non-definite version of matrix orthogonal 
polynomials. They are described by a Riemann-Hilbert matrix Y consisting of four blocks 
Yi,i, Yi t 2, 5^2,1 and ^2,2- In this paper, we show that detYi,! (detY2,2) equals the average 
characteristic polynomial (average inverse characteristic polynomial, respectively) over the 
probabilistic ensemble that is associated to the MOP. In this way we generalize classical re- 
sults for orthogonal polynomials, and also some recent results for MOP of type I and type II. 
We then extend our results to arbitrary products and ratios of characteristic polynomials. 
In the latter case an important role is played by a matrix-valued version of the Christoffel- 
Darboux kernel. Our proofs use determinantal identities involving Schur complements, and 
adaptations of the classical results by Heine, Christoffel and Uvarov. 

Keywords: Multiple/matrix orthogonal polynomials, Christoffel-Darboux kernel, Riemann- 
Hilbert problem, determinantal point process, average characteristic polynomial, Schur com- 
plement, (block) Hankel determinant. 

1 Introduction and statement of results 
1.1 Random matrix ensembles 

On the space Ti. n of Hermitian n by n matrices consider the random matrix ensemble defined by 
the probability distribution 

1 e -TW(M) rfMj M G (L1) 

Z n 

for some given polynomial V of even degree. Here Z n is a normalization constant, Tr denotes 
the trace and dM is the Lebesgue measure on 7i n - 

The random matrix ensemble (II. ip leads to a probability distribution on the space 

R< ={(x 1 ,...,%)el"| xx<...<x n } 

of ordered eigenvalue tuples. It is well-known that this joint probability distribution has the form 



" i<3 J=l J=l 



IP 1 II'' (1.2) 



for some normalization constant Z n . Thus the probability that the ordered eigenvalues of the 
matrix lie in an infinitesimal box [x\ 1 x\ + dx\] x . . . x [x n , x n + dx n ] C M< is given by ljl.2|L 
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Note in particular that the density (|1.2p is small if two eigenvalues xi and Xj are close to each 
other. This means that the eigenvalues tend to 'repel' each other. 
One can write (jl.2p alternatively as 

1 n 

A-det(/ l (x i ))" i=1 det( 5l (x i ))" j=1 J| dx, (1.3) 

where 

/i(i)=a; < - 1 > = x^e"^. (1.4) 

Indeed, this follows upon recognizing (|1.2p to be basically a product of two Vandermonde deter- 
minants. 

Rather than the space R< of ordered eigenvalue tuples, we will find it convenient to con- 
sider the probability distribution (|1.3p - (|1.4j) on the full space K ra . To maintain a probability 
distribution we should then multiply the normalization constant Z n by a factor n\. 

To the probability distribution (|1.3[) - (|1.4p one can associate the average characteristic poly- 
nomial 

-i roc />oo / n \ n 

P n (z) = w \ ... [H(z - Xj ) det (Mxj))^ det Gfcfo))^ U dxj. (1.5) 

J — oo J-oo \j—i I j — i 

It follows from a classical calculation of Heine, see e.g. that P n can be characterized as 

the nth monic orthogonal polynomial with respect to the weig ht function e~ v ^ on M. Thus the 
polynomials P n satisfy the conditions 

P n (x) =x n +0{x n ^) 

for all n 6 N and 

(%')l~ > m (.xjg ^ ^ dx — c n c m J m . n (1-6) 

for all n, m £ N, for certain c n G M. 

Intuitively, the above result states that the zeros of the monic orthogonal polynomial P n (z) 
determine the 'typical' eigenvalue configuration of the random matrix ensemble (jl.ip . This holds 
in particular in the Gaussian case V(x) = x 2 /2, where the ensemble (ll.ip reduces to the well- 
known Gaussian unitary ensemble (GUE) while the corresponding orthogonal polynomials are 
(up to scaling factors) the classical Hermite polynomials. 

In the literature, many more results on average characteristic polynomials for random matrix 
ensembles can be found, see e.g. [11 [3 HI EH [25] and the references therein. We mention the 
following result of Fyodorov-Strahov [16j . Define the average inverse characteristic polynomial 
corresponding to the probability distribution (11.3D - (|1.4p as 

-i ^.oo /.oo / n \ n 

Qn(z) = -j- / • • • / \l[(z - x,)- 1 det (fi(xj))l j=1 det (ftfo))^ \{ dx j} (1.7) 

J-oo J — oo \j = \ I j = i 

for z e C \ K. Then it holds that 

Q n (z) = -±- / ^ n ' l[X)e dx, (1.8) 
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with c„ as in (|1.6|) . Fyodorov-Strahov 16, 25J also observed that both P n (z) and the right hand 
side of (|1.8p have a natural interpretation in terms of the associated Riemann-Hilbert problem 
(briefly RH problem), and they obtained determinantal formulae for averages of more general 
products and ratios of characteristic polynomials; see further. 

In recent years there has been interest in the following generalization of p.l[) , 

1 e -Tr(V{M)-AM) dM _ (LQ) 

Here A is a fixed diagonal matrix which is called the external source. Typically A has only a 
small number of distinct eigenvalues 01, . . . ,a p , say with corresponding multiplicities n±, . . . ,n p . 

The model (|1.9[) was first studied by Brezin-Hikami EH] and P. Zinn- Justin [SU] who 
showed that the eigenvalue correlations are determinantal. In [5] it was shown that the eigen- 
values x\, . . . , x n have a joint probability distribution of the form (|1.3|l . where now fi{x) — x 1 ^ 1 
(Vandcrmonde factor) while gi(x) are given by several Vandermonde-like series: 

9i (x) = ^-i e -Wx)-a 1 x) ) i = l i ... : n 1: 

g ni+i (x) = x i-l e -(V(x)-a 2 x) } j = 1, . . . >n2> 

a , , — T i ~ 1 e~ i - v( - xS, ~ apxS > i — 1 n 

yn 1 + ...+n p - 1 +t \^ J ■*< e , i — x, . . . , lip. 

To this eigenvalue ensemble one can associate the average characteristic polynomial P n (z) in 
exactly the same way as before, see (ll.5[) . Bleher-Kuijlaars [S] showed that P n {z) satisfies the 
multiple orthogonality relations 

/oo 
P n {x)x i e- {v{x) ~ akX) dx = 0, 
-oo 

for i — 0, . . . ,nk — 1 and k = I,..., p. For more background on this kind of orthogonality 
relations, see e.g. 0[2H]. Desrosiers-Forrester [T^] showed that the result ([1.80 on the average 
inverse characteristic polynomial (|1.7p can also be generalized to the random matrix ensemble 
with external source (|1.9| . 

The goal of this paper is to generalize the above results to the more general context of multiple 
orthogonal polynomial ensembles (briefly MOP ensembles) in the sense of Daems-Kuijlaars [llj . 
We will show that in general, (jl.5p and (|1.7p can be expressed as Riemann-Hilbert minors, i.e., 
determinants of certain submatrices of the Riemann-Hilbert matrix. Using these results, we will 
then obtain determinantal formulas for averages of arbitrary products and ratios of characteristic 
polynomials, by adapting the method of Baik-Deift-Strahov [1]. Our results will be stated in 
Section [LU In Sections 1 1 . 2H1 .31 we first recall the basic definitions concerning MOP ensembles. 

1.2 MOP ensembles 

We consider a stochastic model in the following way [TT]. Let p, q £ N be two positive integers. 
Let there be given 

• A (finite) sequence of positive integers n\,ni, . . . ,n v £ N; 

• A sequence of weight functions w\^\(x), w\^,{x), . . . , Wi tP (x) : R — > M; 

• A sequence of positive integers mi , m-i , . . . , m q £ N; 
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• A sequence of weight functions W2,i(x), iV2,2(x), . . . , W2, q {x) : K — ► M. 

We will use the vector notations n := [n\, . . . , n p ), |n| := 2fc=i n & an< ^ similarly for m and |m|. 
Occasionally we will also write wi(x) := (wi t i(x), . . . ,wx >p (x)) and similarly for v?2(x). Assume 
in what follows that |n| = |m| =: n. 

To the above data one can associate a stochastic model, called a multiple orthogonal polyno- 
mial ensemble or briefly MOP ensemble. This model consists of n random points x\, . . . ,x n on 
the real line whose joint p.d.f. can be written as a product of two determinants: 

±r det {fi{ Xj ))l j=1 ■ det fate))" , (1.10) 

where Z n is a normalization factor, and with functions 

fi(x) =x % ~ 1 wx tl (x), i=l,...,m, 
fm+i(x) = x^wi^ix), i = 1, 



,n 2 , 



fm+...+n p -t+i( x ) = x% 1 Wi, p {x), i = l,...,n p , 

and 

gi(x) = x i ~ 1 w 2 ,i(x), i = l,...,mi, 

g mi+i (x) = x t ~ 1 w 2 ,2(x), i = l,.. . ,m 2 , 

9mi + ...+m,_i+iW = X*~ W2,q(x), 1 = 1,..., TUq. 



(1.11) 



(1.12) 



Note that this is a special case of a biorthogonal ensemble [S]. The motivation for the name 'MOP 
ensemble' comes from the multiple orthogonal polynomials that are related to it, see Section [2l 

In order for the above model to be a valid probability distribution on W 1 one should have 
that (jl.lOp is positive on R". We will not address this topic here since our algebraic results will 
be valid irrespective of this positivity condition. 

The normalization constant Z n in (11.10|) serves to make the total probability of the MOP 
ensemble on 1" equal to 1. We will further obtain an expression for Z n as a block Hankel 
determinant; see the remark at the end of Section EOI 

Of course, the above model contains the eigenvalue ensembles in Section [TTT1 as special cases. 
Further motivation comes from the theory of non-intersecting one-dimensional Brownian paths 
with p distinct starting positions and q distinct ending positions, see Figure [TJ See [TT] [20] for 
details. Another type of application can be found in [14) . 

To facilitate comparison with the literature, we may note the following terminology which is 
often used: the MOP ensemble is said to be of type I if q = 1, of type II if p = 1 and of mixed 
type if p, q > 2. Thus the results of our paper will be valid for general MOP ensembles of mixed 
type, whereas the results of [3 [THl H2] only apply to MOP ensembles of type I or type II. 



1.3 Riemann-Hilbert problem and kernel 

The MOP ensemble in Section 11.21 and its associated multiple orthogonal polynomials are in- 
timately related to the following Riemann-Hilbert problem (RH problem) introduced in [llj . 
The RH problem generalizes the well-known RH problem for orthogonal polynomials due to 
Fokas-Its-Kitaev [T3] as well as its generalization in [25] . 
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Figure 1: The figure shows n = 20 non- intersecting Brownian motions with p = 2 starting points 
and q = 3 ending points. The horizontal axis denotes the time t £ [0, 1] and the vertical axis 
denotes the space variable x. For each t £ (0, 1) the positions x\, . . . ,x n of the paths at time t 
form a MOP ensemble (fTTT0l) - (fTTT2]l . 



RH problem 1.1. Consider the data n, Wi(x), m and W2(i) as above and assume that |n| = 
|m| =: n. The RH problem consists in finding a matrix-valued function Y(z) = l^ n ,wi,m,w 2 ( z ) °f 
size p + q by p + q such that 



(1) Y(z) is analytic in C \ M; 

(2) For x £ K, it holds that 



Y + {x)=Y.{x)[^ 



I P W(x) 



In 



where Ik denotes the identity matrix of size k; where 

/w 1 s(x)w 2 ,i(x) ... w ll i(x)w 2 , q (x)\ (w 1 ^i(xY 

W(x) = 



(1.13) 



(W2,l(x) ... W 2 ,q{x)) , 



\W 1 .p(x)w 2 ,l{x) ... Wi )P (x)w 2 .q{x) J \wi^ p (x) / 

(1.14) 

and where the notation Y + {x), Y— (x) denotes the limit ofY(z) with z approaching x G R 
from the upper or lower half plane in C, respectively; 

(3) As z —y oo, we have that 

Y{z) = (I p+q + 0(l/z)) diag^" 1 , . . . , z"-, z-™ 1 , . . . , (1.15) 

The solution Y(z) to the RH problem is unique, if it exists. Partition this matrix as 



p q 
Y(z)= p (Y hl (z) Y h2 (z)\ 
q \Y 2A (z) Y 2a {z) ' 



(1.16) 



where the partition is such that Yi.i(z) has size p x p, and so on. Then the entries of Yi.i and 
Y 2t \ can be described in terms of multiple orthogonal polynomials, while Y12 and ^2,2 contain 
certain Cauchy transforms thereof; see Section [5] for details. 
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The RH problem can be used to define the Christoffel-Darboux kernel. There are actually 
two kernels into play. The first is a q x p matrix- valued kernel K n (x,y) defined as 



K n (x,y) 



2m(y — x) 



1 




(1.17) 



for x, y € C, x ^ y. This definition may seem rather unmotivated; but see Section [3 for the 
connection with the classical definition of the Christoffel-Darboux kernel. 
Second, there is also a scalar- valued kernel K n (x,y) defined by 



It is known that the MOP ensemble in Section 11.21 is determinantal with correlation kernel 
K„(x,y) pi]- In this 

paper, however, the quantity of interest will be the matrix-valued kernel 

K n (x,y). 

1.4 Statement of results 

Now we are ready to state our main results. Throughout this section we assume fixed data n, 
wi(x), m and W2(x) as before, and we consider the corresponding MOP ensemble (|1 . 10p — (|1 . 12|) . 
For any K, L e N we define Pjl ^ clS 



i.e., Pk M is the average with respect to the MOP ensemble (|1.10p - (|1.12j) of the ratio of products 
of characteristic polynomials. Here we assume that y\, . . . , yx £ C, z\,...,zl G C \ M and the 
numbers in the set (yi, . . . , yn, %i, ■ ■ ■ , Zl) are pairwise distinct. 

Note that for K = 1 and L = 0, (|1.19p reduces to the average characteristic polynomial (| 1 . 5(1 . 
while for K = and L = 1 it reduces to the average inverse characteristic polynomial (|1.7jl . 
These are the cases under consideration in our first two main theorems. 

Theorem 1.2. (Average characteristic polynomial) We have that 




(1.18) 



Pl K ' L ^y 1 ,...,y K ;z 1 ,...,z L ) : 




dxi . . . dx n , (1-19) 



P^(y) = detY 1A ( y ), 



(1.20) 



according to the partition (|1.16[) . 



Theorem 1.3. (Average inverse characteristic polynomial) We have that 



P^(z) = det Y 2 , 2 (z), 



(1.21) 



according to the partition Ijl.l6|) . 

The proofs of Theorems 11.2111.31 will be given in Section 13.21 



Remark 1.4. 



G 



(a) One may understand Theorem 11.21 as follows. The quantity detY± t i[y) is a monic poly- 
nomial of degree n, and its zeros determine the 'typical' point configuration of the MOP 
ensemble in Section 1X751 By specializing this to the MOP ensembles in the introduction, 
one obtains the typical eigenvalue configuration of random matrices with external source, 
the typical positions of non-intersecting Brownian motions, etc. However, one should be 
careful with these interpretations. For one thing, we do not even now if all the zeros of 
det Yi,i(2/) are real. This is an open problem. 

(b) The functions P\^ ,L ^ in (|1.19|) are clearly invariant under the involution (n, wi) <-> (m, W2). 
The similar fact for the right hand sides of (|1.20[) and p.21[) is discussed in Section [2J see 

As we mentioned before, Theorem l 1 . 21 generalizes a classical result for orthogonal polynomials, 
see e.g. [12], as well as its generalization to type II MOP due to Bleher-Kuijlaars [5]. On the other 
hand, Theorem l 1 .31 generalizes a result for orthogonal polynomials due to Fyodorov-Strahov [IB] , 
and also its generalization to type I MOP due to Desrosiers- Forrester [T3] , see also [52] . 

Now we state our third main theorem. 

Theorem 1.5. (Average ratio of characteristic polynomials) For the case where K = L = 1 in 
(|1.19j) . we have that 

Pt 1] {y,z)= det L n (y,z) (1.22) 

where L n is the q by q matrix 

L n (y,z) := I q -(z- y) J K n (y,x)W(x)/(z - x) dx = (0 I q ) Y- 1 (y)Y(z) . (1.23) 
Alternatively, we also have 

Pt 1] (y,z)= det Rn(z,y) (1.24) 

where R n is the p by p matrix 

R n (z,y) := I p -(z-y) J°° W(x)K n (x,y)/(z - x) dx = (l p 0)Y- 1 (z)Y(y)( I A. (1.25) 
Here W(x) is the weight matrix (|1.14[) and K n (x,y) is the Christoffel-Darboux kernel (|1.17[) . 



The proof of Theorem 11.51 will be given in Section 13.11 

Note that the matrices W(x) and K n (x, y) are of size p x q and qx p, respectively, so each of 
the matrix expressions in (|1.23|) and (|1.25|) is well-defined. (The integrals are taken entry- wise.) 
The equalities between these expressions are discussed in Section! 



Corollary 1.6. Assume that q = 1 and W2,i(x) = 1. Then 

/•OO 

Pi 1 ^(y,z) = l-(z-y) / K n (y,x)/(z-x) dx, (1.26) 



where K n is the scalar kernel in (|1.18ll . In particular, 

K n (y, z) = — lim -L (P r pl (y, z + ei) - (y, z - ei)) , zel. (1.27) 
z — y e->0 2m V / 



y 

Similar results hold when p — 1 and wi.i(x) 
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Formula (|1.27p retrieves a result in [13] , see also [7] for the scalar case p = q = 1 . 
Finally, let us consider the case of arbitrary K and L in (|1.19p . To this end we need some 
extra notation. 

Definition 1.7. Assume that |n| = |m|. We set (no, mo) := (n, m) and we define the sequence 
of multi- indices (rife, m^), each n& being a vector of length p and each a vector of length q, 
for k G Z inductively as follows 

• For k — 1,2,..., we set n^ := n^-i + (1, 1, .... 1) and we fix arbitrarily such that 
mfc > nifc_i componentwise and — m/j-i) = p. 

• For k = —1, —2, . . ., we set := m^+i — (1, 1, . . . , 1) and we fix n*, arbitrarily such that 
Hfc < nfc + i componentwise and jn^+i — n^ = q. 

Definition 1 1 . 71 implies that | rife | = |rrifc| for each k, so we can consider the RH problem with 
respect to the pair of multi- indices (nfc,mfc). We will assume that the RH problem is solvable 
for all involved pairs of multi- indices. 

Note that in Definition 1 1 . 71 the multi-indices (n^, m^) for negative k are only meaningful pro- 
vided \k\ < mhii—i^ ^q mi. Otherwise we would obtain multi-indices with negative components. 

Now we are ready to consider (|1.19|) for arbitrary K and L. It turns out that the results 
can be expressed as determinants of large matrices constructed from the 'building blocks' in 
Theorems 1 1 . 2flT31 To get the idea we first describe the three typical cases. 

Theorem 1.8. For the expression Pf M in (fTTT9j) . we have 

(a) IfL = then 

n ,m (vk) 



P ^° ] ' • • • ' = n Iv~^¥ det 

\.\.i<i<i<K\yj yi> 



\(Yi,i) (Vi) ■■■ (Yi.i) n K _i,mi(_i 

(1.28) 



(b) IfK = then 

( (^2,2) n , m n ,m \ Z L) 



P^(z u ...,z L ) = = det 



Here we assume that L < 1 + min^ =1 to/ 
(c) If K = L then 

P} l K ' K ^y 1 ,...,y K ,z 1 ,...,z K ) 



\(Y2,2)n_ {L _ 1) , m _ {L _ 1) {zi) ■■■ {Y2,2)n_ {L _ 1) , m _ {L _ 1) {zL) J 

(1.29) 



det 



\Ul<i<j<K(yj ~ Vi)( Z i ~ Z j) J ' 1 s 1 p ', s 

C 

where R n is the matrix in ()1.25p . A similar result holds with the matrix L n in p.23j) 
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Theorem 11.81 generalizes results for the scalar case in [TB] . Our proof will be based on the 
methods of 0]. 

We note that Parts (a)-(b) of Theorem 11.81 are of a different flavor than Part (c). In fact, 
it is shown in [4j [16] that in the scalar case, one can obtain 'mixtures' of the formulas in Parts 
(a)-(b); but these formulas do not seem to have an analogue in the present setting. This is the 
reason why we need to work with the formula of 'two-point' type in Part (c). 

The case of arbitrary K, L in (|1.19p is obtained from mixtures of either Parts (a) and (c) , or 
Parts (b) and (c) of Theorem 1 1.81 



Theorem 1.9. For the expression P, 
(a) If K > L then 

P^ L \ yi ,...,y K ,zx,...,z L ) 
I 



K,L] 



()1.19p . we have 



{-i) L(K - L) Y\a 



det 



Y[l<i<j<K(Vj Vi) Y\l<i<j<h( Z i Z j) / 



-Rn(z L ,yi) 



zl-vi 
0l,l)n ,m (2/l) 



J^Rn{z L ,y K ) 
(ll.Ono.mofe) 



(1.31) 



Wx,x) (yi) ■■■ <Xi,i) {Vk)J 

Here all R n matrices are taken with respect to the pair of multi-indices no, mo 
(b) If L> K then 



I 



det 



T\l<i<j<K(yi Vj) Yil<i<j<L( Z j Z i) ) 



j^L n (y K , Zl ) 

(Zl) 



j^L n (y K ,z L ) 

n , m 



V(^2, 2 )„_ ( 



L-K-l),m-(L-K 



-l)( Z l) 



(^2,2)n_ (i __ K _ 1 



(1.32) 



) ,m-(L-K-l) 



Ml 



Here all L n matrices are taken with respect to the pair of multi-indices no, mo, and we 
assume that L — K < 1 + min^ =1 to/. 

1.5 About the proofs of Theorems 11.51 and 11.8( c) 

Theorems 11.51 and II. 8f c) will be the key results from which all other main theorems will follow. 
The proofs of these two theorems will be based on a Schur complement formula for the kernel 
K n (x,y). This formula is valid for an arbitrary weight matrix 



W(x) = (^,iW) fc=1 pM 



(1.33) 
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of size p by q. Thus W(x) does not need to have the rank-one factorization (II . 14|) . For any 
k e {1, . . . ,p} and I G {1, . . . ,q} let 

/OO 
-OO 

denote the j'th moment of the weight function Wk,i{x), and for any r, s € N let 

rr(M) _ f„(WA 

V ' • /j=0,..,r-lj=0,.,.-l 

be the r x s Hankel matrix formed out of these moments. Stack these matrices in the block 
Hankel matrix 

/fr(l-l) rr(l,9) \ 



H n 



Note that this matrix is of size n x n. 
Recall that for a block matrix 



l rr(P,l) 



(1.34) 



rr(P>9) I 



M := 



A B 
C D 



(1.35) 



with A square of size k x k (say), the Schur complement of M with respect to the submatrix D 
is defined as the matrix 

S m ,d := D — CA- X B, (1.36) 

provided A is invertible. Schur complements are also known as quasi- determinants in the litera- 
ture [18j . For more information on Schur complements see Section ll.61 

Proposition 1.10. The kernel K n (x,y) equals the Schur complement of the matrix 



' ^ni.mi 



(1,9) 



til ,1Tlg 



l± rip ,m ( 



V 



o \ 



(1.37) 



/ 



with respect to its bottom right q x p submatrix. Here the superscript T denotes the matrix 



transpose, and we use the column vector notations y„ 



1 y ... y m 1 ) and x. 



(1 x 



Proposition ! 1 . 1 01 will be proved in Section l275l Variants of this proposition in different contexts 
can be found in [6l [23] , see also [3] . 

Proposition II . 101 will allow the quantities det L n (y, z) and deti?„(z, y) in Theorem 11.51 to be 
written as a ratio of two determinants. In Section T3. II we will expand these determinants using an 
adaptation of the classical argument of Heine, and this will lead us to the proof of Theorem 1 1.51 
In Section [3.31 we will give a similar argument proving Theorem II .8( c). 
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1.6 Background on Schur complements 

Since Schur complements frequently occur in this paper, we find it convenient to list here some 
basic properties. 

Recall the setting in (| 1 . 35|) and (|1.36p . Schur complements are related to Gaussian elimination 
as follows: 

'A B\ (A 



C Dl - V S A - ■ I (L38) 



where the ~ symbol relates matrices that can be obtained from one another by multiplying on 
the left and on the right with square transformation matrices of the form ( k z\ and 



* I J \0 I j 

respectively. The procedure to move from the left to the right hand side of (|1.38p is sometimes 
called Gaussian elimination with pivot block A, see |19j . 

In the special case where D and (hence) Sm,d are square matrices, one obtains from (|1.38|) 
the determinant relation 

detS M , D = ^-. (1.39) 
det A 

Let Si j denote the (i, j)th entry of S = Sm,d- By ()1.36p one has 



</,, C.l : b,. (1.40) 



where c, denotes the ith row of C and hj denotes the jth column of B. Now (|1.40|) is nothing 
but the Schur complement of the matrix 

A b 3 



(1.41) 

with respect to the entry dij. From (|1.39p we then obtain the determinantal expression 

det 



A bj 



det A 



(1.42) 



This shows that the entries of the Schur complement are ratios of determinants. 

More generally than p.4ip , one may observe that Schur complements are well-behaved with 
respect to taking submatrices, in the sense that any submatrix of the Schur complement (|1.36p 
is itself a Schur complement, of an appropriate submatrix of (jl.35p . 

Finally, observe that for any matrix U of appropriate size, the matrix product U Sm,d — 
U(D — CA~ 1 B) can again be written as a Schur complement, of the matrix 

A B \ M AnS 

UC UDj- (L43) 

A similar fact holds for matrix products of the form Sm,dV. 



1.7 Outline of the rest of the paper 

The rest of this paper is organized as follows. In Section [2] we discuss auxiliary results and 
notations that are used in the proofs of the main theorems, and we prove Proposition 11.101 
Theorems 11.51 and ll.8f c) are then proved in Section by an adaptation of the classical argu- 
ment of Heine. Theorems 11.21 and 11.31 are obtained as limiting cases of Theorem 11.51 Finally, 
the generalizations to arbitrary products and ratios of characteristic polynomials are proved in 
Section [H 



11 



2 Preliminaries for the proofs 



In this section we collect some auxiliary results and notations that are used in the proofs of 
the main theorems. This section is organized as follows. Vector orthogonal polynomials (which 
include multiple orthogonal polynomials as a special case) are defined in Section \2. II and their re- 
lation to the RH problem in Section l2~2l The connection with block Hankel determinants and the 
duality relations are discussed in Sections 12.31 and 12.41 In Section 12.51 we discuss the Christoffel- 
Darboux kernel, leading to the proof of Proposition ll.101 Finally, the related quantities L n and 
R n are investigated in Section [2.61 

Remark 2.1. Practical note: the reader who is mainly interested in Theorems 11.21 [T31 and 1 1 . 51 and 
who wants to take Proposition 11.101 for granted may skip this entire section and move directly 
to Section [31 



2.1 Vector orthogonal polynomials 

Let p, q G N and consider the multi-indices n, m as in the beginning of Section 11.21 but now 
with |n| = |m| + 1. Let W(x) be a p by q weight matrix as in (|1.33|) . Our definition of vector 
orthogonal polynomials will be the same as in Sorokin-Van Iseghem [24] ; this definition includes 
multiple orthogonal polynomials as a special case, when the weight matrix W(x) has the rank-one 
factorization (|1.14|) . 

Let P n be the following space of polynomial vectors 



r n = {P( x ) = 



with Pfc(x) a polynomial of degree at most rik — !}• 



(2.1) 



The standard basis (Ai(a;), . . . , A| n |(x)) of the vector space V n is defined by the column vectors 
of the p by |n| matrix 



1 x 




(Ai(a:) ... A N (a;)) = 

\ 

Similarly, let V m denote the vector space 





1 x 







r ri2— 1 



\ 





(2.2) 



V m = {Q(x) = 



with Qi(x) a polynomial of degree at most mi — 1}, 



(2.3) 



\Qq( X )/ 



and define the standard basis (Bi(x), . . . ,B| m i(x)) for V m by the columns of the q by |m| matrix 



(Bi(a;) ... B w (a:)) = 



/ 1 x 




V 



„mi — 1 





1 x 







r m 2 — 1 



\ 





(2.4) 



Note that we use boldface notation to denote vector-valued objects. 
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We say that P(x) = P n ,m(^) G Tn is a vector orthogonal polynomial with respect to the 
multi-indices n, m and the weight matrix W(x) if 

/CO 
Pl m (x)W(x)Q(x)dx = 
-co 

for every Q e V m . To stress the dependence on the weight matrix we will sometimes write 

Pn,m(s) = Pn,m,w(x). 

The coefficients of the vector orthogonal polynomials can be found from a homogeneous linear 
system with |n| unknowns (polynomial coefficients) and |m| equations (orthogonality conditions). 
The restriction |n| = |m| + 1 guarantees that this system has a nontrivial solution. If P nim (a;) 
is unique up to a multiplicative factor then the pair of multi-indices n, m is called normal. 

Assume that n, m is a normal pair of indices. Then P ni m(x) is said to satisfy 

• the normalization of type I with respect to the Ith index, I £ {1, . . . , q}, if 




Pl m (x)W(x)Q(x) dx = 1 



where Q(x) is the vector consisting of zeros except for the Ith. entry which is x mi . 

• the normalization of type II with respect to the kth index, k 6 {1, . . . ,p}, if the fcth com- 
ponent (Pk)n,m(x) is monic, i.e., if 

The vector orthogonal polynomials P n .m(a;) corresponding to the above normalizations, if 
they exist, will be denoted as Pn,m(x) and Pnfm\x), respectively. 



2.2 Solution to the RH problem 

Consider again the RH problem in Section [T~3l with |n| = |m| but now with W{x) in (| 1 . 14[) 
replaced by an arbitrary matrix (|1.33p . The solution to the RH problem, if it exists, is uniquely 
described by vector orthogonal polynomials and their Cauchy transforms [llj . More precisely, 
consider again the partition of Y(z) as in (|1.16|) . Then one has that 

• Y\ y i{z) has its fcth row given by the row vector 

k = l,...,p. 

• ^2,1 (z) has its lih row given by 



-2™(P^_ ei ) T (z), (2.6) 



l = l,...,q. 

Yi^(z) has its fcth row given by 



2tti 

, ,p. 



(.II, k) \T ( 

dx, (2.7) 
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• ^2.2 (z) has its Ith row given by 



(pi!lj T ww(x) 



dx, 



(2.8) 



1 = 1, 



Here we are using the notation e; to denote the standard basis vector which is zero except for 
its Zth entry which is one. The length of this vector will always be clear from the context. 

For example, in the special case where p = q = 2, the solution Y(z) is given by the 4x4 
matrix 

/ CP 



Y(z) = D x 



(-Pi) n+e2 ,m 
CP 

V 1 /n,m— ei 



(P 

v 2/n+ei . in 
(p\(H,2j 
\ 2 /n+e2 ■ in 

V ^ / n . in - e i 

.m — 02 



where D := diag(l, 1, 
as in (l2~7l) and K 



(P 2 p 

-2m, —2iri), and where the entries denoted with * are Cauchy transforms 



2.3 Biorthogonality and moment matrix 

Recall the standard bases Ai(x) and Bj(x) in ()2.2|1 and (12. 4|) . where now |n| — |m| =: n. The 
moment matrix M with respect to these bases is of size n by n and has entries 



Mij = 



/OO 
Af(x)W{x)B j {x) dx. 
-oo 



The moment matrix coincides with the block Hankel matrix H n m in ()1.34p . 

Consider ordered bases (Pi(x), . . . , P„(i)) and (Qi(aj), . . . , Q„(x)) of the vector spaces V n 
and V m , respectively. We say that these bases are biorthogonal if 

/oc 
Pf(x)W(x)Q j (x) dx = S id , (2.9) 
-oo 

the Kronecker delta. 

Lemma 2.2. Assume that [n| = |m|. Then the following statements are equivalent 

1. The RH problem for Y{z) in Section ] is solvable. 

2. There is no non-zero vector in Vn which is biorthogonal to the entire space V m . 

3. dct# n , m 5* 0. 

4- There exist biorthogonal bases for the spaces V n and V m . 

PROOF. First consider the equivalence between Statements 1 and 2. Suppose that statement 2 
holds true. Then each of the polynomial vectors in l|2.5p - (|2.6p exists [TT] and therefore the RH 
problem for Y(z) is solvable. Conversely, suppose that statement 2 does not hold. Then the 
polynomial vectors in (|2.5|) - (|2.6|) . if they exist, are not unique and therefore the RH problem for 
Y{z) cannot be solvable, since that would contradict the uniqueness of Y{z). Finally, the equiva- 
lence between Statements 2-4 follows from standard linear algebra arguments whose description 
we omit. □ 
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2.4 Duality 

The role of the biorthogonal bases Pi and Qj is interchanged by swapping the multi-indices n 
and m and by transposing the weight matrix W(x): 

(n, m, W(x)) <-> (m, n, W T (x)). (2.10) 

This duality has a nice form in terms of the RH problem [TJ [TT] . If we partition the RH matrices 
corresponding to the original and dual data in (|2.10p as 

PI 3 P 

Y n , m , w = P (Y\,i Y ia \ Y w t = q (Y lfX Yx <2 



q \Y2,i Y2.2) 1 p \ Yi,\ Y 2 .i 

respectively, then it holds that 



Y\a Y\2\ / Y2.2 —Y21 



Y 2 ,i Y2.2) \-Yia *i, 



(2.11) 



where the superscript ~ T denotes the inverse transpose. From a well-known theorem on minors 
of the inverse matrix, see e.g. [TTl page. 21, eq. (33)], and the general fact that dety(z) = 1 it 
then follows that 

dety 1 , 1 (7/) = detF 1 , 1 ( 2 /) (2.12) 

and 

detY 2>2 (z) = dety 2 , 2 (z). (2.13) 

The duality relations (|2.12() ()2.13|) take care of a point made earlier, see the remark after the 
statement of Theorem 11.31 

From another application of [T71 loc. cit.] it also follows that 



det ((I p 0) Y- 1 (z)Y(y) = det ((0 I,) Y^{y)Y{z) ^ 



(2.14) 



Comparing this with (|1.23[) and (| 1 . 25[) . this establishes the fact that det R n (z, y) — det L n (y, z), 
which is of course implicit in the statement of Theorem 11.51 



2.5 The Christoffel-Darboux kernel 

In this section we prove the Schur complement formula for the kernel K n (x, y) in Proposition ll.ldl 
To this end we identify K n (x,y) as a reproducing kernel. In this way we will show that (|1.17p 
corresponds to the kernels familiar from the theory of matrix orthogonal polynomials |23j and 
vector orthogonal polynomials [24] . 

We will assume throughout that the RH problem for Y(z) is solvable, which is of course 
necessary in order for (|1.17p to make sense. 

Lemma 2.3. The kernel K n (x,y) is a q x p matrix whose entry is a bivariate polynomial 

of degree at most mi — 1 in x and rij — 1 in y. Equivalently, 

71 

K n (x,y) = J2 OijBi(x)Aj{y), (2.15) 
for certain Cjj, where Aj and denote the basis vectors in (|2.2p and (|2.4p . 
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PROOF. The expression (0 7 g ) Y 1 (x)Y(y) \^qJ m the numerator of p,17|) is made out of 

the polynomial entries of the matrices Y~ 1 (x) and Y(y), cf. (|2.5[) - (|2.6[) and (|2.1ip . Since this 
expression vanishes if x = y, each of these polynomials is divisible by y — x, so the denominator 
of (|1.17|) can be divided out. □ 

It was shown in pj] that for a rank-one weight matrix W(x) as in (|1.14[) . the scalar- valued 
kernel K n (x, y) in (|1.18[) has a certain reproducing property. We now obtain a similar property 
for the matrix-valued kernel K n {x,y) (without restrictions on W(x)). 

Proposition 2.4. The kernel K n satisfies the reproducing property 

/oo 
K n (x, y)W(y)Q(y) dy = Q(x) (2.16) 
-oo 

for any i£C and Q 6 P m . Moreover, this reproducing property uniquely characterizes K n (x,y) 
in the class of bivariate polynomial matrices of the form (|2.15|) . 

PROOF. First we establish (|2.16jl . This can be done by adapting the proof in [11]. However, 
we give a more streamlined proof. We have 



K n (x,y)W(y)Q(y) dy = J K n (x,y)W(y)(Q(y) - Q(x)) dy + / K n (x,y)W(y)Q(x) dy. 
By inserting (|1.17[) this becomes 



MJY 2 Ay)w(y)^E^ 



= (0 ^-^TfY^^dy)^ vy-^&x^Z)™ 

(2-17) 

The first term in (|2.17|) is zero since j s a polynomial vector in y whose Zth entry has 

degree at most m; — 2, I — 1, . . . , q, and by invoking the orthogonality relations of the vector 
orthogonal polynomials. For the second term in (|2.17p . we recognize the defining relations for 
the Cauchy transforms in the second block column of the RH problem, see (|2.7p - (|2.8p . so we get 



Q(aO = (0 I q ) ( r ° ) Q(x) = Q(x). 



This establishes (|2.16p when x G C \ M; the case where i6l follows by continuity. 

Next we show that the reproducing property (|2.16[) uniquely characterizes K n (x,y) in the 
class of bivariate polynomial matrices (|2.15p . By Lemma [2~2l 1 =>• 4, we can choose biorthogonal 
bases for the polynomial vector spaces V n and V m ; let (Pfe)^ =1 and (Qfc)fc=i be such bases. Any 
q by p matrix M n (x, y) of the form (|2.15|) can be rewritten as 

n 

M n {x,y) = c it jQi(x)Pj{y) 

for suitable constants Cjj. But then the reproducing property (|2.16p and the biorthogonality 
relations (|2.9p imply that Cjj = Sij, the Kronecker delta. This ends the proof. □ 



Proposition 12.41 has the following dual version. 
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Proposition 2.5. The kernel K n satisfies the reproducing property 

P T (x)W(x)K n (x,y) dx = P T (y) 



(2.18) 

for any y £ C and P £ "P n . Moreover, this reproducing property uniquely characterizes K n (x,y) 
in the class of bivariate polynomial matrices of the form (|2.15[) . 

Now we are ready for the 

PROOF of Proposition 11.1(71 Denote with M n (x,y) the Schur complement of Q1.37|) . We 
check the reproducing kernel property 



A T (x)W{x)M n (x,y) dx = A T (y), 



(2.19) 



for each column vector of the form A(x) = (0 ... x lk ... 0) with i^ £ {0, . . . , — 1}, 
k £ {1, . . . ,p}. (These vectors A(x) are the standard basis of the space V n -) By using the 
multiplication property (|1.43[) for Schur complements, the left hand side of (|2.19p equals the 
Schur complement of 



Hi 



(1,9) 



Til ,Triq 



H, 



(p,i) 



H 



i h (*,i) 



(M 
ik,m q 



o \ 



o J 



(2.20) 



with respect to its bottom right 1 xp submatrix, where we use the row vector notation r .. 
( f^ k + - ) • Subtracting from the last row of (|2.20|) the ifcth row from the fcth block row, 

\ %k 3 / j=0,...,mi— 1 

we are led to the Schur complement of 
/ ttW 



rr(p,l) 
1J -n p ,mi 



Tj(p,q) 



\ 



o 



o 







-v 



(2.21) 



/ 



where now the last row is zero except for the entry — y %k in the fcth column of the second 
block column. But clearly, the Schur complement of (|2.21[) is just the row vector A T (y) = 
(0 ... y lk ... 0) . This establishes the reproducing kernel property (|2.19[) . 

To finish the proof of the proposition, we note that the degree structure of (|1.37[) implies 
that M n (x,y) has the form in (|2.15p . Thus the proof is ended by invoking the uniqueness part 
of Proposition l2~5l □ 



2.6 The matrices L n and R n 

In this section we consider in more detail the quantities L n and R n which are derived from 
the kernel K n . Let us first establish the equivalence between the two different formulae in the 
definition of L n in (|1.23p . From the identity y — x = (y — z) + [z — x) we obtain 

^—^K n {y,x)W{x)dx = [ K n (y,x)W(x)I q dx + [ ^^K n (y, x)W(x) dx 
z — X J J z — X 

= I q -(z-y) [ K n (y,x)W(x)/{z-x) dx (2.22) 
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by virtue of the reproducing property in Proposition ^. 41 On the other hand, the equality 

/ JZ^K»(v,x)W(x) dx = (0 I q ) Y- 1 (y)Y(z) (?) 

follows from (fl~T7| and (j^TT ]) — ([2T5 ]l . This establishes the required equality in (TPS]) . 

The equivalence between the different formulae in (|1.25p can be established similarly. 
In section SI we will need the following property of L n . 

Proposition 2.6. The matrix L n in (|1.23p satisfies the 'vanishing property' 

f 00 V T {y)W{y)^^dy = Q (2.23) 

7-oo z-y 

for any z£C\I and P G V n . 

PROOF. From the expression of L n in (j2.22[) . we obtain 

[p-(y)W(y)^l d y= dy - / / * T {v)W{v)K n {y, X )W{ X ) ^ 

J ' z-y J z-y J J z-x 

which is 

= J p T (y) EMl dy-J(J P T (y)W(y)K n (y,x)dy^ ^dx. 
By inserting the reproducing property this becomes 

= [pT(y)^dy-[rT {x) ^ dx = 0. 
J z-y J z-x 

□ 



Here is the dual version of Proposition 12.61 
Proposition 2.7. The matrix R n in (| 1 .25(1 satisfies the 'vanishing property' 



r ^^y} W{y) Q {y) dy = (2.24) 

J-oo z-y 



for any z £ C \ K and Q6?„ 



3 Proofs, part 1 

In this section we prove Theorems II. 2\ [L3l 11.51 and ll.8f c). 
3.1 Proof of Theorem [131 



First we prove Theorem 11.51 The proof of this theorem will follow by expanding the moment 
determinant in Proposition II . 101 in a similar way as in the classical argument of Heine. 

We will restrict ourselves to the proof of (|1.24|k the proof of p.22|l can then be devised 
in a similar way, or by simply invoking (|2.14|) . For notational convenience and to keep things 
readable, we give the proof for the case p = q = 2. The case of general p, q is completely similar 
except that it requires more notational burden. 
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First we expand the left hand side of (|1.24[) . By definition, the average ratio of characteristic 
polynomials Pn' X \y, z) is given by 



3=1 



v~ x 3 



W\,l{Xl) 



x" 1 ... xH 1 1 w 1 , 1 (x n ) 



X T 1 wi. 2 (xi) 



X„ 2 1 Wio(Xr, 



W2,l(xi) 



W2,l(x n ) 



x'i 1 w 2 ,i{xi) ... x™ 1 1 w 2 ,i(x n ) 



W2,2i x l) ■■■ W 2 ,2(x n ) 

x" l2 ~ 1 W 2 ,2(xi) ... X™ 2 ~ 1 W2,2{Xn) 



\\ dx 3 

3=1 



(3.1) 



We can expand the second determinant according to the Lagrange expansion with respect to the 
first block row: 



w 2 ,i{xi) 



W2,i[x n ) 



XT 1 ~ 1 W 2 ,1(X 1 ) ... C 1 " 1 ^,!^) 



W 2 ,2(xi) 



W 2 ,2(Xn) 



xT 2 1 W 2 ,2(xi) ... X™ 2 1 W 2 ,2{x n ) 

W2,l(x Sl ) ... W 2 ,l{x Smi ) 

X^- 1 W 2 .l{x Sl ) ... X™l- l W 2 ,l{x Sm J 



W2.2{x Sl ) 



W2,2(Xg ) 



717, o — 1 / \ Wlo — 1 / \ 

X H W2.2{Xs 1 ) ... X s 2 W2,2(X Sm ) 



(3.2) 

where the sum is over all (£ ) subsets S C {1, ...,n} with \S\ — mi, where we write the 
elements of S in increasing order: si < . . . < s mi , similarly for those of the complement S — 
{1, . . . , n} \ S, Si < ... < s m2 , and where we denote by (— l) s the sign of the permutation 

For each term of ()3.2() we can incorporate the factor (— l) s into the other determinant of (|3.ip 
by permuting its columns. Next we can relabel the variables Xi,...,x n back in their original 
form, and so (|3.1[) reduces to 



W2,l(xi) 



TT V ~ x 3 

*-\ Z -Xj 

3 = 1 J 



W2,l(x mi ) 



wi,i(x n ) 



-^"V.ifri) ... x^- 1 w 1 ^{x n ) 



Wl,2(xi) 



Wl,2(x„) 



xT'^lMxi) •■■ xl 2 ~ l W la {x n ) 

W2,2(x mi + l) ■■■ W 2 ,2(x n ) 



X? 1 W2,l(xt) ... xZ\~ 1 W2,l{x mi ) 



xZ\+\w2,2(x mi + 1 ) ... X™ 2 1 W 2 ,2(Xn) 



]Jdx 3 

3=1 
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which on evaluating the Vandermonde determinants becomes 



if' 11 = 



(") 

\mi/ 



n 



wi,i(x„) 



x" 1 wi,i(xi) ... x^ 1 1 w 1 , 1 (x n ) 



Wl,2(si) 

712 — 1 / \ 

x x 101,2(2:1) 



ioi, 2 (x n ) 
x™ 2_1 u;i i2 (x„) 



Kt<j<mi mi+l<i<j<n J — 1 j-mi + 1 j — 1 



Next we expand the right hand side of (|1.24|) . We recall from Proposition ll.lOl that — K n (x, y) 
is the Schur complement of the matrix 



(3.4) 



17(1.1) 

n 1 ,mi 

17(2,1) 

Eini ,mi 


77(1,2) 
77(2,2) 


y«i 




T 
X mi 






T 
X m 2 







/ 



This implies that jz^^2 — J w ( x ^ n J x ' v * > dx is the Schur complement of the matrix 



z— y 
ff (i,i) 

1J -ni ,mi 

u-(2,i) 

-^712 ,mi 



o-(l,2) 

1J -n\ .7712 

77(2,2) 

- n «2,™2 



J^i^liW/l 2 - J" x m 2 Wi,2(2:)/(2 - x) da; 

V I^m^iixyiz-x) dx Jx^ 2 W 2 ,2{x)/(z-x) dx 



From p.39p this implies the determinant representation 









y« 2 



1/(2-1/) 

l/{z-y) J 



(3.5) 



det 



Z - y 



W(x)if»(x,y) 



rfx 



(1,1) 

Tii ,mi 
(2,1) 



77(1,2) 

11 1 .777,2 

17(2,2) 

- rz Ti2,m 2 







JxmjWuW/^- 1 )^ /x^ 2 Wi, 2 (a;)/(z-x) dx l/(z-y) 
Jx^ 1 W 2jl (x)/(z-x) dx Jx^ 2 T^ 2 , 2 (x)/(z-x) dx 





y« 2 
o 

1/(^-2/) 



where k = deti/ n!m is the determinant of the top left block of (|3.5|) . 

Now we expand the integral in the right hand side of (|3.6|) . By definition, each H 



(3.6) 



r (k,0 



is a 



Hankel matrix whose («, j)th entry is the (i + j)th moment with respect to the weight function 

Wkj(x) = wi > k(x)w 2 ,i(x), for k,l = 1,2: 



i) v . 



x i+j W k ,i(x) dx. 



(3.7) 



Note that the integration variable x in (|3.7p is just a 'dummy variable' which we can give any 
name we want. By choosing the same integration variable xi for all entries in the first column 
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of (|3.6p , we can rewrite (|3.6p as follows: 





fWi ;1 (xi) dx x 


* 




* 


* 


1 







J xT^W^xx) dx x 




* 


* 


* 


ynx-x 





1 


J W 2 ,i{x 1 ) dx 1 


* 




* 


* 





1 


K 


J a^Wb.lfrl) dX! 


* 




* 







y n 2 -l 




f Wl - l(Xl) d Xl 
f W2 - l{xi) dxx 

J Z — X\ 1 






* 
* 


* 
* 


1 

z-y 




1 

z-y 



where the entries denoted with * are not important at this stage of the proof. Using the multi- 
linearity of the determinant with respect to its first column, we can take out the integration with 
respect to x%: 



Wx,l(xx) 


* 


* 


* 


* 


1 





a;" 1_1 wi,i(a;i) 


* 


* 






y m-i 





wi,2(a;i) 


* 











1 


< 2 -V,2(zi) 


* 


* 









y n 2 -l 


z—x\ 


* 


* 






l 

z—y 





Z — X\ 


* 


* 









1 

z-y 



W2A(xi)dx± 



Here we also used that W\,x{x) = wi i i(x)u>2.i( x ) and W2,i(x) = W\ t 2{x)w2,i(x), allowing us to 
extract the common factor W2,i(xi) from the first column. Applying the same technique to the 
other columns of the first block column we obtain 



wx,i{xx) 


Wx,l(Xmx) 


* 


* 


1 





X^^Wi^Xi) . 




* 




yUX-l 





Wx,2(xx) 


Wx,2(Xmi) 


* 


* 





1 


XT~ X W1,2{XX) . 




* 


* 





y n 2 -l 


Z — Xi 


z — x mi 


* 


* 


1 

z-y 





Z — X\ 


nil, i(x mi ) 
z — x mi 


* 







1 

z-y 



Y\ (x 3 f 1 w 2 ,i(x j )d2 



3 = 1 



To get rid of the factor rij=i x f ' one can a PPly the following well-known symmetrization trick. 
Sum the above expression over all permutations x ai , . . . , x„ m of Xx, . . . , x mi . For each term in 
this sum, apply a column permutation to arrange the columns of the determinant back in their 
original form; this releases a factor (— l) a . But then we get the expansion of a Vandermonde 
determinant: 

mi 

EH) ff n< i= n fe-o- 

<J 7 — 1 Ki<j<mi 
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So we get 



/■c(mi)! 



101,1 (xi) 


u>i,i(x mi ) 




* 


1 













y ni-l 
















1 













y n 2 -l 


Wl,l(xi) 






* 


1 





Z—X\ 


z—x mi 






*—v 


l»l,2(ll) 


«Jl,2(x mi ) 









1 


Z—X\ 


z— x mi 
















rj 1 1 


i 



l<i<j<mi 

Applying the same techniques to the second block column, we find 



n ~ x ^> n ^A^) n 

3 = 1 3=1 



K(mi)!(m 2 )! 



ioi,i(a;i) 




1 





x r l 1 ~ 1 w lt i(xi) . 




y m-i 





wi t2 {xi) 


wi t2 (x n ) 





1 


x 7 l 2 ~ 1 wi t2 (xi) . 


. x^ l 2 ~ 1 wi t2 {x n ) 





y n 2 -l 


WX,l(Xx) 


wi : i{x n ) 


i 





Z — X\ 

uii,a(a:i) 

2-^1 


z—x n 
z—x n 


*—v 



1 

*-v 



n ( x 3~ x i) n ( x o - x i)Y[ w 2,i(xj) n ^2,2(^0 jj^. (3.8) 

j=l i=mi+l j=l 



l<i<_j'<mi mi+l<i<j<n 

We claim that this can be rewritten as 



1 



1 



(z — y) 2 K,{mi)\(m 2 )\ 



n^ifa-z.?) 

TV=i(z-Xj) 



wi,i(x n ) 



"T_ 1 wi^(x 1 ) ... x^ 1 1 w 1 , 1 (x n ) 



wi i2 (xi) 



wi i2 (x n ) 



x[ 2 1 wi^ 2 (xi) ... x™ 2 1 wi, 2 {x n ) 

mi n n 

n ( x 3~ x i) n ( x 3 - x i)Y[ w 2,i( x j) n ^a^yi^- ( 3 - 9 ) 

l<i<j<mi m\+l<i<j<n j — 1 j— tni+1 J — 1 

This follows from standard determinant manipulations and the exact formula for a Cauchy- 
Vandermonde determinant; see the appendix. 

Now up to a constant factor times l/(z — y) 2 , (|3.9|) equals (|3.3| . This proves (|1.24[> . up to a 
constant factor; but this constant factor must be 1 as can be seen by matching the leading terms 
on both sides of (|Qgj) . 
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Remark 3.1. By matching the leading constants in (|3.3p and (|3.9p . we obtain the following 
expression for the normalization constant Z n in (jl.lOp : 

Z„ = nnl 

where we recall that k is the block Hankel determinant k = det H n m . This identity also follows 
directly from a 'generalized Cauchy-Binet identity', see e.g. [55]. 



3.2 Proof of Theorems Q and IT51 

Having established Theorem 11.51 Theorems II .21 and II. 31 now follow as simple limiting cases. For 
Theorem 11.21 we let z —> oo in (|1.24[) and observe from (|1.25p that 



detiU^y) = det ((l p 0) Y~ 1 (z)Y(y) (^Jj = det ^((Y^) 1A (z) (Y~ 1 ) 1:2 (z)) & 



YiAv) 
(y) 

z- n det Y 1A (y) + 0(z- n - 1 ), 



where in the last step we used the Cauchy-Binet identity and the asymptotics (| 1 . 1 5[) . This 
establishes Theorem 11.21 Theorem 11.31 can be obtained in a similar way by letting y — > oo in 
CL22J). 



3.3 Proof of Theorem fTTM c) 

The proof that we have given for Theorem 1 1.51 can be extended to prove Theorem I 1.8f c) as well. 
The key observation is that the matrix in the right hand side of (|1.30[) can be written in Schur 
complement form as well. For example, for K = L = 2 and p = q = 2 one has that 



is the Schur complement of the matrix. 



#(1,1) 

11 n± ,mi 

rr(2,l) 

- n -7i2,rrti 




rr(l,2) 
J - L m -rii2 
rr(2,2) 






(yi) 




rii 




(yi)na 


(y 2 ) 




"1 


\ 

(y2)« 2 


/ W M (*)/(*! 


— x) dx 


fxl 2 W h2 (x)/( Zl 


-x) 


dx 


i/( Zl - 


-yi) 







V(*i - 


- V2) 







/x^Wa,! (x)/(zi 


— x) dx 


/x^Wa, a (x)/(«i 


-x) 


dx 







V(^i - 









1/0*1 - 


2/2) 


J^ ni W 1A (x)/(z2 


— x) dx 


/x^Wi, a (x)/(«2 


-x) 


dx 


1/0*2 - 


-wi) 









- 3/2 ) 







f^ ni W 2A (x)/(z2 


— x) dx 


S* T m2 W 2a {x)l(z 2 


-x) 


dx 







1/0*2 - 









l/(*2 - 


2/2) / 



(3.10) 



Compare with (|3.5|) . 

The determinant of this Schur complement can be manipulated in the same way as in Sec- 
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tion l3.ll This leads to the following analogue of (|3 .8[) : 



1 



K(mi)!(m2)! 



wi,i(xi) 



Wi t i(x n ) 



x" 1 ... x™ 1 1 wi A (x n ) 



Wi i2 (xi) 







1 





1 


no — l / \ 

x t 2 Wia(Xl) . 


. X^^WiaiXn) 





no — 1 

Vi 





no — 1 

y 2 



t"l,l(ftl)' 



21- 


11 




2=1) 


zi- 




i«l,l 


2=1) 


Z2- 


Kl 


U)l,2 


2=1) 


^2- 


2=1 



n ~ x ^ n ( ; 



Z\ — x n 
zi—x n 

Z 2 —X„ 

mi 



Hi — 1 n ni — 1 



2\— yi 










l 

^2— yi 



Zl-J/2 


1 

Z2-J/2 





21— 2/2 


1 

Z2-J/2 



l<i<j<mi mi+l<j<j<n 

and the following analogue of (|3.9|) : 



(21 - z 2 )(y 2 ~ Vi) 



Xi)X\_W2,l{Xj) Y[ w zA x ]) H dx ji (3-11) 
j=l j=i 



nL=ife-%) 7 k(™i)K™ 2 )< 



l<z<j<mi mi+Ki<j<n 



iui,i(xi) ... wi,i(x„) 
^ 1_1 ioi,i(xi) ... £c" 1 ~ 1 u;i > i(a; n ) 



Wl, 2 (xi) 



X 



/ / / 1 



1 1 



n 



Xi) JJ W2,i(xj) I I lOa.afe) JJetej. (3.12) 

j=l j=mi+l j=l 



Here the transition from (|3.11jl to (|3.12| follows again from standard determinant manipulations 
and the exact formula for a Cauchy-Vandermonde determinant, see the appendix. 

Formula (|3.12[) gives an expression for the determinant in the right hand side of (| 1 . 30[) . Here 
the leftmost factor 

/ {zi - z 2 )(y 2 - Vi) \ 



Ili,,-=i(*i ~Vi) 



(3.13) 



is canceled by multiplying with the first factor in the right hand side in (|1.30p . Finally, by also 
expanding the left hand side of (|1.30|) in the same way as before, cf. (|3 . 3|) . Theorem ll.8f c) follows. 



4 Proofs, part 2 

In this section we prove Theorems ll.8f a)-(b) and 11.91 about averages of general products and 
ratios of characteristic polynomials. To this end we use a similar strategy as in Baik-Deift- 
Strahov [3], i.e., we make use of the already established Theorems 11.2 1 [T~5l together with adapta- 
tions of the classical results by Christoffcl and Uvarov. We mention that these Christoffcl-Uvarov 
type results will be valid for arbitrary vector orthogonal polynomials, i.e., the weight matrix W(x) 
does not need to have the rank-one factorization p,14[) . 
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Throughout this section we assume that |n| = |m| and we let W(x) be an arbitrary weight 
matrix of size p x q. We let Y(z) be the solution of the corresponding RH problem and K n (x, y) 
the kernel (|1.17|) . We also consider the modified weight matrix 

W(x) = ^= liX ~ Vk) W(x), (4.1) 
Ui=i( x ~ z i) 

and we denote the corresponding RH matrix and kernel with Y(z) and K n (x, y). Here we assume 
that . . . , yx € C, Z\, . . . , Zl 6 C \ R and the numbers in the set (y\, . . . , yx, Zi,..., z£) are 
pairwise distinct. Our aim will be to relate the matrices Y and Y, and similarly for the kernels 
K n and K n . 

Although Theorem 1 1.81 is just a special case of Theorem 1 1.9[ for reasons of comprehension we 
discuss the former first, see Sections 14.11 and 14.21 The proof of Theorem 11.91 is then only briefly 
sketched in Section |4~31 

4.1 Proof of Theorem Of a) 

Recall that Christoffel 's formula expresses the orthogonal polynomials with respect to a polyno- 
mially modified weight function w(x) = rife=i( a: — Vk)w{x) in terms of those with respect to the 
original weight function w(x) [H[5S]. We now s ^ 3 ^ e a generalization of this formula to the context 
of vector orthogonal polynomials. Throughout, wc will use the notations in Definition 11.71 

Proposition 4.1. Assume that L — in (|4.1[) . so that 

K 

W(x) = Y[(x-y k )W{x). (4.2) 

fe=i 



Then Y\^(y) is the Schur complement of the matrix 



1 


( Ol,l)n ,nio(j/l) 


(Xl,l)n 




(>l,l)no,mo(y) 




(yi,l)n K -l,m K -i(yi) • 


•• (^l,l)n K -l 


m e -i(») 


(Y 1 ,i) riK _ 1 ^ mK _ 1 (y) 




V (Y l .i) YlK . mK {y l ) 


{Y\,l)a K 


m K (yi<) 


(Yl,l)n K ,m K {y) 



(4.3) 

with respect to its bottom right p by p submatrix. 

PROOF. Denote with S(y) the Schur complement of the matrix (|4.3|) . ft is clear that S(y) is 
a matrix of size p by p. We claim that the entries of S(y) are polynomials in y, and moreover 
that 

S(y) = (l p + 0(l/y))d^g(y n \...,y^), y -> oo. (4.4) 

To prove these claims, first note that all entries in the last block column of (|4.3[) are polynomials 
in y. Now we use the expression (jl.42l) for the entries of the Schur complement of the matrix 
(ll.35[) . It is clear that if we apply this formula to the input matrix (|4.3p , then the determinant 
in the numerator of (jl.42|) is zero whenever y — yu for certain k = 1, . . . , K . This means that 
all entries of S(y) are polynomials divisible by Jlfc=i(2/ — 2/fe)i which on multiplying with the 
prefactor 1/Ilfc=i(y — Dk) in (|4. 3[) indeed shows that S(y) is a polynomial matrix. 

To establish (14. 4p . observe that the leading term of S(y) comes from the bottom right p by 
p block in (|4.3|) . which has the asymptotics 

(Fi.iW ,m K (y) = (I P + 0(l/y)) diag(y"i +K , . . . , y n *+ K ), y -> oo. 
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Thus on dividing by the prefactor 1/ IlfcLi^ — Vk) m H-^P - we obtain (|4.4|) . 
Next, we check the orthogonality conditions 



S(x)W(x)Q(x) dx = 



(4.5) 



for all Q 6 V m , where we recall the definition of W(x) in (|4.2[) . To prove this, note that the 
factor n^i( a; — Vk) m 114. 2|) cancels with the prefactor in (|4.3j) . and by the multi-linearity of 
determinants the remaining integral can be taken into the last block column of (|4.3|) . So the left 
hand side of (|4.5[) is the Schur complement of the matrix 



/ (Yi.i) 



1,1 Jno,mo 



(1/0 



0T,l)n o ,m o (2/if) 



f(Y 1<1 ) Il0iIao (x)W(x)Q(x) dx \ 
J(Y hl ) nK _ umK _ 1 (x)W(x)Cl(x) dx 



n K _i,m K _i 

j(Y 1 , 1 ) nK . mK (x)W(x)Q{x) dx J 

(4.6) 

But each of the matrices Y\ t i in the last column of ()4.6j) satisfies the orthogonality relations with 
respect to Q 6 P m . So the last column of (|4.6j) is zero and therefore the Schur complement of 
this matrix is zero as well. This establishes (|4.5I) . 

The proof of the proposition now follows from (|4.4p and (|4.5p . see (|2.5[) . □ 



Having established Theorem 11.21 and Proposition 14. 1[ and taking into account the property 
(jl.39|) . Theorem 1 1.8( a) can now be proved by induction on K, in exactly the same way as in [4j 
Proof of Th. 4.3]. We refer to the latter reference for details. 

4.2 Proof of Theorem fTMb) 

Recall that Uvarov's formula expresses the orthogonal polynomials with respect to an inverse- 
polynomially modified weight function w(x) — w(x) / Yif=i ( x ~ z i) m terms of those with respect 
to the original weight function w(x) [U[27]. We now state a generalization of this formula to the 
context of vector orthogonal polynomials. 

Proposition 4.2. Assume that K = in (|4.ip . so that 

W(x) = — — W(x). 



T\i=iix ~ Zl) 

Then la.iO^) is the Schur complement of the matrix 
( 



(4.7) 



(^2,2)no,mo( z l) 


(>2,2)n 


m„(zL) 


0"2,l)n 


m ( z ) 


(^2,2)n_ I . + i,m_ I , + 1 (2 : l) • 


• (^2,2)n_z, + i 


m__L+i ( z l) 


(^2,l)n_ £+1 


m_ i + i ( z ) 


(5 A 2,2)n_t,m_ z ,(^l) 


(^2,2)n_ £ 






m_x,0) 



(4.8) 



with respect to its bottom right q by p submatrix, and ^2,2(2) is the Schur complement of the 
matrix 



rL=iO - z i) 



(^2,2)n ,m (2 : l) 


(^2,2)no 




(^ / 2,2)n ,m ( 2; ) 


(^2,2)n_ I/+ i,m_ I , + i(zi) . 


• • (^2,2)n_ i + i 


m_ i+ i { z h) 


(l2,2)n_ Il+ i 1 m.i + i(«) 


(^ / 2,2)n_ £ ,m_ i (^l) 


(^2,2)n_i 




(i2,2)n_x,,m_ I , (*) 



(4.9) 



26 



with respect to its bottom right q by q submatrix. 

PROOF. Denote with S(z) the Schur complement of the matrix (|4.8|) . It is clear that S(z) is 
a polynomial matrix of size q by p, which has the required degree structure 

S(z) = (0(l/z))diag(z ni ,...,z n "), z^oo. (4.10) 

On the other hand, denote with T(z) the Schur complement of the matrix (|4.9|) . From the 
asymptotic condition 

(^2,2)n_^m_ t (z) = (I q + 0(1/ z)) diag^^™ 1 , . . . , Z^), 

it is clear that T(z) has the required asymptotics 

T(z) = (I q + 0(l/z))diag{z- m \...,z- m «), z^oo. (4.11) 
Taking into account (|4.10|) and l|4.11[) . the proposition will follow if we can show that 

i r s(x)w(x) 



T{z) 



2ni 



dx. 



Consider the left hand side of (14.12|) . From the formula 12,2(2) = ™ / 



Y 2 ,i(x)W(x) 



(4.12) 



dx we see 



that all entries of matrix (|4.9|) are integrals. Choose a common integration variable xi for all 
entries in the «th row. Consider the coefficients in the partial fraction decomposition 



L 

E 

1=1 



Q 



1 



X - Zl 



{x-z)n l=1 (x-zi)' 



(4.13) 



where q, c are functions of {z\ 1 . . . ,zl, z} but not of x. Consider the column operation where to 
the last block column of matrix (|4.9p we add ci times the Ith block column, I = 1, . . . , L. Perform- 
ing these operations inside the integrals, this causes the integrand in the last block column to be 
multiplied with c/ Yit=i( x ~ z i) = ~ z i)l T\lLi( x ~ z i) ■ Here the factor n^i( z ~ z i) can be 

taken out of the integrand, canceling the prefactor in (|4.9p . while the factor W(x)/Ylf—i(x — zi) 
in the integrand is nothing but W(x). So we obtain precisely the right hand side of (|4.12p . This 
establishes (|4.12[) . and thereby the proposition is proved. □ 



Having established Theorem 11.31 and Proposition 14. 2\ and taking into account the property 
(ll.39p . Theorem 1 1.8T b) can now be proved by induction on L, in exactly the same way as in [H 
Proof of Th. 4.10]. We refer to the latter reference for details. 



4.3 Proof of Theorem PTT91 

Finally we prove Theorem II. 91 The key to the proof of Theorem [TT^a) is the following result. 
Proposition 4.3. Assume that K > L in (|4.ip . Then Yl,! (j/) is the Schur complement of the 



matrix 



( 



nf=iO/ - z i) 



-Rn(zi,yi) 



ir ^ r i?„(z L ,y 1 ) 
(>l,l)no,mo (2/l) 



V (Yl,l) 



n K - L ,m K - L 



(yi) 



-R n {z L ,y K ) 



zl-Vk 
(il,l) no ,m (VK) 



l,l;njf_i,mK- 



-Rn{zx,y) \ 



-R n (z L ,y) 



zl-v 

(Yi,i) no , mo (y) 



(Vk) (Yi,i) 



n K -L,m K - L 



(y) I 

(4.14) 
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with respect to its bottom right p by p submatrix. Here all R n matrices are taken with respect to 
the pair of multi-indices no, nio. 

The proof follows the same lines as the proof of Proposition 14.11 and we omit the details. 
We suffice to mention that the proof makes use of the vanishing property for the matrix R n in 
Proposition 12.71 

From Proposition I4.3[ Theorem 11.9( a) can now be proved by induction on K, K > L (for L 
fixed), with induction basis K = L corresponding to Theorem ll.Sf c). 

Finally, we note that the proof of Theorem 1 1.9f b) hinges on the following fact. 

Proposition 4.4. Assume that L > K in (|4. 1[) . Then Y 2 ,i (z) is the Schur complement of the 
matrix 



n , m 



J^LniVuZL) 



(zl) 



n K - L ,m K 



-Azl) 



2-KiK n (yx,z) \ 
2niK n (y K , z) 

(^2,l)n ,m (^) 



(4.15) 



(z) J 



with respect to its bottom right q by p submatrix, and 12,2(2) is the Schur complement of the 
matrix 



Uk=i( z -yk) 

Uti( z - z i) 



7^ L n{VK,Zl) 
n , m 



j^L n {y K ,z L ) 

(^2,2)no,m (^) 



V (Y 2a ) 



njc-l,,m 



k-l( z i) 



(zl) 



^L n (y u z) \ 



j^L n (y K ,z) 



(^2,2)n ,mo(^) 



(Y 2 ,2)n K - L ,Tn K _ L (z) J 

(4.16) 

with respect to its bottom right q by q submatrix. Here all L n matrices are taken with respect to 
the pair of multi-indices no, niQ. 



PROOF. The proof follows the same lines as the proof of Proposition 14.21 The properties 
(|4.10|) and (|4.1ip are easily extended to the present situation. The main difficulty lies in proving 
(|4.12p . To this end one can use the following analogue of the partial fraction decomposition 

6131), 



E 



Q 



1 



(4.17) 

i=i - - - - z )Ui=i( x - z i) 

where q, c are functions of {yi, . . . , yic, Zi, . . . , Zl, z} but not of x. Observe that these very same 
coefficients ci , c also appear in the more complicated partial fraction decomposition 



E 



ci(x - y k ) 



x - yk 



cllfcLl^ - Vk) 



(4.18) 



f-[ (zi - yk){x - zi) (z-y k )(x-z) (x-z)Uf =1 (x-ziy 

for any k. The proof of (I4.12|) then follows as in the proof of Proposition ^. 21 using the integral 
representations Y 2i2 (z) — — ^ / ^^^^ dx (as before) and 

L n (y,z) = / - — -K n (y,x)W(x) dx, 
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and applying the relations (14. 17[) for the integrals for ¥2,2 and (|4.18|) for those for L n . □ 



Appendix: Calculations with Cauchy-Vandermonde deter- 
minants 



Define the functions 



fiix) = 



i = 1, . . . , n, 



l/{zi- n -x), i = n+ 1, . . . ,n + m, 
where z\, ■ ■ ■ , z m are a given set of numbers. It is well-known that 

n+rn _ Yll<i<j<m( Zi ~ Z i) Y[l<i<j<n+m( X j ~ X i) 



det(/,fe)):;™ 



(4.19) 



The expression in the left hand side of (|4.19[) is called a Cauchy- Vandermonde determinant. 
Proof that $£E§ implies < K9\t 

By applying a suitable row permutation, the determinant inside (|3.8|) takes the form 



(-1) T 





wx,\{x n ) 


1 





z—x-i 


■ x™ 1 ~ 1 w 1 . 1 (x n ) 

z — x n 


1 

«-» 








wia(x n ) 





1 


Xi 2 ~ 1 w 1 . 2 {xi) . 

u)i, 2 (ai) 

Z — X\ 


■ x^ 2 ^ 1 w 1 . 2 (x n ) 

t»l,2(l») 

z — x n 






y n 2 -l 

1 

z—y 



By applying Lagrange expansion with respect to the first block row, this equals 
wi.i(x Sl ) ... wi tl {x Sni ) 1 



= £(-i) f 



X™ 1 ^1,1 

t"i,i(^ si ) 



<- n ^ 1 w 1 , 1 (x Sni ) y 

wi,i(x Sni ) 
z—x Sni 

wi t2 (x Sl ) 



m — 1 



1 



t»i, 2 (a;5 1 ) t"i,2(^a„ 2 ) 1 

z-x Sl • ■ ■ «-i in , z—y 



(4.20) 



where the sum is over all ) subsets 5 C {1, . . . , n} with IS 1 ) = m, where we write the elements 
of S" in increasing order: si < . . . < s ni , similarly for those of the complement S = {1, . . . , n} \S, 
Sx < ■ ■ ■ < s, l2 , and where we denote by (— l) s the sign of the permutation si , . . . , s ni , si , . . . , s n2 . 
Note that the factor (— l)™ 2 of the previous formula has disappeared. 
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Up to a diagonal scaling of the columns, both determinants in (|4.20|) are Cauchy-Vandermonde 
determinants. By virtue of (|4. 19[) . this yields 



n 



y~ x 3 



wi,i{x Sl ) 
x n s l- l w 1A (x Sl ) 



(4.21) 



— 1 / \ no— 1 / \ 

But this is precisely what one gets from Lagrange expansion of the determinant in (|3.9|) . 
Proof that fLTTTD implies fETT2D 



The proof that (|3.11| implies (|3.12|) follows in the same way. The only difference is that there 
is an extra row and column in each Cauchy-Vandermonde determinant, thereby yielding the 
pre-factor (|3. 13[) . 
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